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SUMMARY 

The interactions of a vortex wake with a turbulent strati- 
fied atmosphere are investigated with the computer code "WAKE." 
It is shown that atmospheric shear, turbulence, and stratifica- 
tion can provide the dominant mechanisms by which vortex wakes 
decay. 

Computations have included the interaction of a vortex wake 
with a viscous ground plane. Here, the observed phenomenon of 
vortex bounce is explained in terms of secondary vorticity pro- 
duced on the ground. This vorticity is swept off the ground 
and advected about the vortex pair, thereby altering the classic 
hyperbolic trajectory. When compared with the decay of a vortex 
pair out of ground plane interaction, it is shown that strength 
of the vortex pair is rapidly reduced by the scrubbing on the 
ground . 

The phenomenon of the solitary vortex is explained as an 
interaction of a vortex with crosswind shear. Here, the vortex 
having the sign opposite that of the sign of the vorticity in 
the shear is dispersed by a convective instability. This in- 
stability results in the rapid production of turbulence which in 
turn disperses the smoke marking the vortex. 

It is shown that the effect of stable stratification on the 
descent of a vortex pair is to halt the descent of the pair. 
Numerical experiments have allowed a simple relation to be de- 
rived which predicts the maximum descent of a vortex pair in a 
linearly stratified atmosphere. 

Lastly, the "WAKE" code has been updated, including new 
algorithms to improve accuracy and speed. In particular, the 
original primitive variable formulation has been replaced by a 
stream function-vorticity formulation, thereby eliminating pre- 
vious numerical problems with regard to maintaining divergence- 
free flow. 



INTRODUCTION 


In Ref. 1 (hereafter denoted I) a two-dimensional, unsteady, 
finite difference, computer code ’’WAKE" was documented and used to 
investigate vortex interactions and turbulent transport in aircraft 
wakes. The thrust of I was to investigate the phenomenon of merg- 
ing or coalescing of vortices in multivortex pair wakes. One major 
conclusion of this study was that vortex pairs can be made to 
interact in a favorable way so as to reduce the hazard induced on 
an encountering aircraft. After this interaction was complete, 
however, it was shown that the decay of the wake in the absence of 
atmospheric effects was a very slow process. The study contained 
herein is an extension of I and was undertaken to determine the 
effects of atmospheric turbulence, shear, and stratification as 
well as to investigate the dynamics of vortex interaction with the 
ground plane. In addition, general improvements were made to the 
"WAKE" code which reduced computational time with no sacrifice in 
accuracy. The study contained herein.; assumes that vortex interac- 
tions are complete, and that vortex merging has occurred and that 
the vortex wake can now be characterized by two counter-rotating 
vortices of circulation T r> with separation b . 

The fact that our atmosphere can be sheared, turbulent, and 
stably stratified simultaneously has made difficult the interpreta- 
tion of vortex data obtained during flight tests. Apparently it is 
quite difficult to measure atmospheric parameters sufficiently to 
quantify the dissipative ability of the atmosphere, particularly 
since these measurements must not necessarily be made near the 
ground. For these reasons, it was decided to analytically investi- 
gate atmospheric effects by devising computational models which 
examined atmospheric turbulence, crosswind shear, and stratification 
separately. It is felt that this approach should both produce a 
maximum understanding as well as quantitative predictions of the 
dissipative ability of the atmosphere. 

This report is organized as follows. In Section 2, changes 
made to the "WAKE" code are detailed including a description as to 
how a viscous ground plane has been augmented into the code. In 
Section 3, results are presented as to how atmospheric effects in a 
neutral atmosphere result in decay of a vortex pair. It is noted 
here that much of the material contained in Section 3 is included as 
an appendix and has been presented at the American Institute of 
Aeronautics and Astronautic ' s 16th Aerospace Sciences Meeting in 
Huntsville, Alabama. In Section 4, vortex descent in a stable atmo- 
sphere is investigated. In Section 5, results of several additional 
computations are presented which include roll up of a vortex sheet 
as well as computations of vortex ascent in an enclosed wind tunnel. 
Lastly, in Section 6, conclusions and recommendations are offered. 
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induced drag 
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Brunt-Vaisala circular frequency, 
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r.m.s. turbulence level 
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circulation 

vortex sheet strength 
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2. THE "WAKE" CODE 


The "WAKE" code is a two-dimensional, unsteady, incompressible 
computer program developed to simulate the dynamics of multiple air- 
craft vortices in a turbulent stratified atmosphere. The two-dimen- 
sional approximation is justified for these flows since gradients 
in the flight x direction are very small compared with changes in 
the vertical z and lateral y directions. Therefore, a numeri- 
cal simulation is equivalent to observing a wake flowfield in a 
plane which is both perpendicular to the aircraft's flight direction 
and fixed to the ground. Incompressibility is justified since the 
Mach number of the swirling velocities in the wake is small. The 
"WAKE" code is operational on the NASA Langley computers . 

2.1 Equations of Motion 

The continuity, momentum, and energy equations govern the solu- 
tion to well-posed incompressible, turbulent, stratified flow prob- 
lems which are of interest in the study of vortex wakes. Unfortu- 
nately, full numerical unsteady solution to these equations are 
beyond today's current computer hardware and software capabilities. 
To overcome this problem, turbulent transport in fluids is estimated 
with transport models which have been developed to be compatible 
with today's resources. The approach used here is to numerically 
solve for ensemble averaged mean field and ensemble averaged second- 
order correlations. The current model is one which has been used to 
predict successfully both quantitative and qualitative features of 
many turbulent flowfields (Refs. 1-4). 

The model equations are briefly reviewed . Let each 
dependent variable be written as the sum of a mean and a fluctuating 
part (e.g., U-^ = U-^ + ui) . After formal substitution into the 
continuity, momentum, and energy equations and ensemble averaging, 
the equations for the mean fields become; 
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We have assumed that the flow is incompressible but have per- 
mitted small changes in the density due to changes in temperature. 
The Boussinesq approximation has been used whereby the only effect 
of the density variation is in the body force term in the mean mo- 
mentum equation. 


In or der to solv e f or the second-order correlation ensembled 
averages, u^u-j and U£0 , which are needed in the mean equations, 
the mean equations, Eq. (2.2) and Eq. (2.3), are subtracted from 
the complete equations. This operation leaves equations for u-^ 
and 0 wh ich ca n b e premultiplied and averaged to yield equations 
for ujU| and u-^0 . These equations are then closed by modeling 
third-order as well as certain second-order correlations. The de- 
tails of this closure are not given here but are: discussed in Ref. 
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where q = u.u. (the summation convention is implied) and A is 
a measure of the macroscale of the turbulence. The constants v c , 
b , A and s are constants of the invariant second-order closure 
model. They are assigned the values 0.3 , 0.125 , 0.75 and 1.8 , 
respectively. These values have been chosen from computer simula- 
tions of turbulent flows where experimental data are available 
(see Ref. 5 for a complete description). The scale A is a mea- 
sure of the large eddies to be found in a flow. The scale is, 
therefore, itself determined from a differential scale equation-? 
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where the constants sq , S 2 » S 3 , and S 4 have the values -0.35 , 
-0.6 , 0.375 , and 0.8 . Equations (2.1) through (2.7) are the 
second-order invariant model of turbulent transport. 


2.2 Superequilibrium 

The above turbulent transport model has a simple limit which 
is useful analytically to fix background turbulent atmospheres 
through which vortices will be allowed to descend. This limit 
exists when second-order correlations have equilibrated to the mean 
field; hence the name superequilibrium is used. The superequilibrium 
limit of the turbulent model is Cfourid by neglecting convection and 
diffusion terms in the equations for the second-order correlations. 
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In addition, the high Reynolds number limit of the dissipation 
terms is taken. The resulting equations for the turbulent correla- 
tions are 
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where it has been assumed that Ui = [U(z),0,0] and that 0 = 0(z) . 
The solution of Eq. (2.8) through Eq. (2.14) may be computed if the 
following definitions are introduced. 
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Substitution of the definitions, Eq. (2.15) and Eq. (2.16) into 
Eq. (2.8) through Eq. (2.14) results in a system of nonlinear alge- 
braic equations which may be solved to give expressions for the 
turbulent correlations in terms of the stratification 30/3z , the 
shear 3U/3z , and the scale A . A discussion of the stabilizing 
effects of stratification and rotation in the superequilibrium 
limit is given in Refs. 6-8. 

2.3 Solution Algorithm 

The solution algorithm is the same as that described in I with 
one significant change; the mean equations have been solved in 
terms of stream function-vorticity variables. Note that Eq. (2.1) 
in the two-dimensional approximation is 


3V , 3W n 
3y 3z 

which allowsthe definition of a stream function ' ip. : 
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a Poisson equation 


- n (2.19) 


Equations' (2.2) with i = 2,3 are then cross differentiated and sub- 
tracted to eliminate pressure to yield a transport equation for the 
streamwise component of vorticity 

re n ‘ ^ - + 

( 2 . 20 ) 

The resulting system of equations finally consists of the 
mean vorticity equation, Eq. (2.20), and Eq. (2.3) through Eq. (2.7). 
The Poisson equation, Eq. (2.19), is then solved for jp using a fast 
direct solver developed by Swartztrauber and Sweet (Ref. 9). The 
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velocity field may then be computed from the definition of the 
stream function. Use of the stream function-vorticity system elimin- 
ates the need to compute the pressure field and automatically in- 
sures that continuity, Eq. (2.17), is satisfied, as opposed to the 
algorithm employed in I where a correction for continuity must be 
added at each time step. In addition, boundary conditions for the 
stream’ function (needed for. solving’ the Poisson equation) may be 
directly related to the velocity components V and W , while the 
choice of boundary conditions for pressure in I is less clear. 


2.4 Boundary Conditions 

The "WAKE" code is designed to accept boundary conditions of 
the form Ae + B(3e/3n) + C - 0 on the variable e on any of the 
four sides of the computational domain. A , B , C may be functions 
of y , z, t and 3/3n is the normal derivative. Two special cases 
of the above are described below. 

Rigid Wall - Law of the Wall 

Vortices descending in ground effect must experience a de- 
intensification as a result of viscous scrubbing along the ground. 

In order to quantify the amount of dissipation which results, a vis- 
cous no slip boundary condition is needed. Instead of resolving the 
wall boundary layer with a finite difference mesh, it is more effi- 
cient to develop an asymptotic wall layer condition. This is accom- 
plished by noting that -asymptotically the velocity shear in the 
neighborhood of a wall (say at z = 0) of surface roughness z q , 
is given by ° 


3v = l 

3z K z 


( 2 . 21 ) 


where U is a to be determined friction velocity and k is von 
Karman's constant. In this wall layer, since (3/3z) >> (3/3y) , 
Eq. (2.21) defines an approximation to the vorticity, i.e., 
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These expressions are obtained by requiring that i|j = ij; 2 = 0 at 
2 = z Q • 

The finite difference grid will not adequately resolve the 
dependent quantities in the steep turbulent boundary layer region. 

We will instead use the expressions, Eq. (2.22) - Eq. (2.24), to de- 
fine a boundary condition at a reference level z r above the wall, 
but within the asymptotic region where these equations are valid. 

The reference level z r is the first mesh point in the computa- 
tional domain. To define boundary conditions, consider the finite 
difference formulation for q and V in terms of ip : 


-ip 


= 


v = 


r+l 


+ 2\p 


(Az) 
~TKz 


'P 


r-1 


^r+1 ‘ *r-l 


(2.25) 

(2.26) 


where Az is the mesh size and ip r represents the streamfunction 
at z r and ^ r +i is the streamfunction at the mesh points above 
and below z r tzr-1 is a fictitious point). Equation (2.26) can 
also be written as 
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Elimination of <P r _i from Eq. (2.25) and Eq. (2.27) yields 


n r = 


„ - 2 <Vn ' 't'r 5 


( 1+ 2 & lo S f ) 


(2.28) 


the boundary condition for n in terms of ip . Eliminating U 
from Eq. (2.22) and Eq. (2.24) yields 
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the boundary condition on ip . 
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Radiation or Advection Conditions 

Vortices descending into a stably stratified fluid will gen- 
erate internal waves which will propagate to the edges of the com- 
putational domain. If upon reaching the boundary, some condition 
is not employed to allow these waves to pass out of the domain, 
they will be reflected back and possibly spoil the computation. No 
one universal technique has been devised to avoid this problem. 
Specification of boundary conditions on a finite computational do- 
main representing an infinitely large stratified domain is both 
challenging and frustrating , requiring both physical insight and art. 
Previously, one technique which we have employed was to apply a 
porous liner near the computational boundaries . This is accomplished 
by adding damping terms to the equations of motion,, the' strength of 
the damping increasihg as the Ijoundhry’ is approached .(Ref. 9) / . This 
technique has been abandoned since it is difficult to accurately 
assess how much artifical damping is actually added to the compu- 
tation. Modest success has been achieved with the use of a 
Sommerfield radiation condition which is described below. 

Boundary conditions of the type 


ft + c H = 0 (2 - 30 > 

propagate the quantify e at velocity c in the n direction. 

If c is the phase speed of internal waves as they approach the 
boundary, Eq. (2.30) should allow these waves to propagate out of the 
domain. The success of this condition rests with the accurate spe- 
cification of c and with a narrow width of the wave number spec- 
trum which must be passed out of the domain. 

If c is a fluid speed, Eq. (2.30) is an advection condition. 
For example, in the case of a pair of descending vortices, a down- 
ward moving coordinate system can be imposed which will keep the 
vortices of interest in the center of the domain. This means that 
fluid will be washing out through the top of the domain. A suitable 
condition for the top boundary would be to set c equal to the 
local vertical velocity when e is a variable which will be advected 
out, e.g., the vorticity. 
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3. VORTEX DECAY IN NEUTRAL ATMOSPHERES 
VORTEX/GROUND PLANE INTERACTION 

The investigation of atmospheric effects as they influence 
the decay of aircraft wake vortices can be separated as to whether 
the atmosphere is neutral or stably stratified. In this section, 
results obtained for descent of aircraft vortices in a neutral 
atmosphere are summarized. Details of the work undertaken may be 
found in the Appendix which is a preprint of the work entitled, 

"The Role of Atmospheric Shear, Turbulence and a Ground Plane on 
the Dissipation of Aircraft Vortex Wakes." This work was presented 
at the 16th Aerospace Sciences Meeting in Huntsville, Alabama. 

Results of three types of computations are presented and discussed. 
They are; 

1. Computations undertaken to evaluate the effectiveness 
of background ambient turbulence on the dissipation of 
a vortex wake. 

2. Computations undertaken to evaluate the effectiveness 
of crosswind shear on the dissipation of a vortex wake. 

3. Computations undertaken to evaluate the effectiveness 
of the ground on the dissipation of a vortex wake. 

The following conclusions have been drawn; 

Vortex pairs in a turbulent background - Numerical computations 
have shown that the decay of a vortex pair in a turbulent atmosphere 
is greatly enhanced by the presence of turbulence. Direct computa- 
tion has shown that the decay of circulation in a square contour 
positioned on the vortex "center" decays as when turbulence is 

present. This is compared to a computed t~l/2 decay when the 
ambient turbulence is zero and the pair diffuses as a result of its 
own internally generated turbulence. Computations have shown that 
turbulence levels typical of the atmosphere can reduce the circula- 
tion by a factor of three or so one minute after passage of a B747 
aircraft. The inability to directly correlate real data with com- 
putations is tied to failure to quantify the scale of turbulence in 
the atmosphere. Work is currently under way to validate models 
which predict integral scale lengths. 

Vortex pairs in a crosswind shear - During flight tests, vor- 
tices are marked with smoke for visualization purposes . Under cer- 
tain conditions the smoke in one vortex would vanish, leaving what 
appeared to be a solitary vortex (Ref. 10). It was conjectured at 
the time of observation that the solitary vortex may pose a more 
serious problem than a pair wake and that the phenomenon was in some 
way brought about by atmospheric shear. Computations run with cross- 
wind shear essentially confirm that shear disperses the vortex whose 
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vorticity is of the opposite sign of the vorticity associated with 
the shear. However, direct computations cannot confirm that the 
remaining vortex will be a more significant hazard than a compar- 
able pair wake. Basically, the decay of the solitary vortex is 
strongly influenced by the turbulence associated with the cross 
shear. It has been speculated, that over tens of minutes it is possible 
to have shear without turbulence and turbulence without shear when 
the atmosphere is stably stratified. Under these conditions, it 
may indeed be possible to have a solitary vortex which may persist 
for a long period of time. 

Vortex ground plane interaction - Vortices in close proximity 
to the ground must feel some decay of strength as a consequence of 
scrubbing of fluid along the ground. To quantify this effect com- 
putations were run with a viscous turbulent law of the wall boun- 
dary condition applied at ground level. The results of these com- 
putations were quite surprising in that the vortices did not follow 
the classical hyperbolic trajectories but rebounded from the ground 
plane. This rebounding or bouncing phenomenon results from secon- 
dary vorticity generated at the ground being swept up and convected 
about the primary vortices. This result is significant in that it 
provides a mechanism by which vorticity of opposite sign can be 
mixed with the vortex wake and thereby augment the decay process . 

Direct computation suggests that rolling moment induced on a fol- 
lower aircraft may be reduced by a factor of 3 after one minute of 
the passage of a B747 aircraft. 



4. STABLE ATMOSPHERIC STRATIFICATION 


In this section the descent of a vortex pair in a stably stra- 
tified, incompressible atmosphere is investigated. This study is 
the first which includes turbulent transport directly. Motivation 
has come from the several previously published analytical investi- 
gations, (Refs. 11-15); most of which are in disagreement. This 
disagreement, as discussed by Widnall (Ref. 16), results from the differ- 
ent accounting of fluid impulse and energy and the entraining/ 
detraining of mass and vorticity from the interface of the descending 
oval of fluid. The computations discussed below will suggest that 
the dynamics of the descent of a vortex pair in a stably stratified 
flow is intimately tied to the correct accounting of turbulent trans- 
port effects. 

411 Physical Effects of Stratification 

Consider a linearly stratified stable incompressible fluid with 
the unperturbed density given by 



(4.1) 


The characteristic scale height H of the stratification is defined 


H = 



(4.2) 


and is assumed to be much greater than the characteristic vortex 
spacing b . The characteristic time of a gravity wave is known 
to be 


x = 
g 



(4.3) 


A measure of the importance of stratification is to compare the 
characteristic vortex pair time x v = 2iTs2/r with the B.V. period. 
This ratio of time scales defines a Froude number which is 
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(4.4) 


The Froude number for a standard atmosphere at an altitude of 
300 m can be computed to be of the order of 10 for a B747 aircraft. 
This implies that stratification effects are weak. An alternate 
interpretation is that the effects of stratification are significant 
after the vortex pair has descended; say five spans. ' The Froude number 
is infinite when the flow is neutrally stratified. 
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It is appropriate here to discuss qualitatively the anticipated 
effects of stratification prior to presenting numerical results. In 
Fig. 4.1(a) is shown the streamline pattern of a vortex pair descend- 
ing in a neutral fluid. As is well known the vortex pair transports 
mass downward. This mass is that contained in the oval of recircu- 
lating fluid which is observed in a coordinate system moving with 
the vortex pair. As is also well known, this descent velocity is 
approximately r/4irs . In the absence of three-dimensional instabi- 
lity and ambient atmospheric turbulence, a vortex pair propagates 
with little change in velocity (see Fig. 2 in the Appendix). If 
this oval of fluid of density p Q were to descend into a fluid of 
density pq such that pq > p 0 °it is anticipated that the behavior 
of the pair is significantly altered as a result of the vorticity 
which will be generated at the density discontinuity. 


The origin of this vorticity can be seen by examining Eq. (2.20) 
where density p may be substituted for temperature 0 , since 
p/p = - 0/0 . The production term is 


i© =_£_ iP 

0 o ay P o 3y 


(4.5) 
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and is proportional to the gradient of the density in the horizontal, 
y , direction. This density gradient is large at the oval interface. 


It should be anticipated that the additional vorticity gener- 
ated will alter the descent of the pair as a consequence of inducing 
an additional component of velocity on the descending vortices. We 
have schematically illustrated in Fig. 4.1 three possibilities. The 
centroids of the buoyancy generated vorticity are above the pair, 
below the pair or at the altitude of the pair. The first-order con- 
sequence of these centroid locations in the absence of turbulent 
diffusion would be to decrease the separation between vortices in 
the pair and possibly lead to an acceleration (Ref. 14) , increase 
the separation or hold constant the separation, respectively. Un- 
fortunately, turbulent diffusion alters this simple description. 


4.2 Descent of a Vortex Pair in a Stably Stratified Atmosphere . 

In order to study the decay of vortex pairs in a stably strati- 
fied atmosphere, simulations were run at three different Froude num- 
bers with three levels of background turbulence (determined from 
superequilibrium theory), nine cases in all. 

The computations are initialized by positioning a Gaussian of 
vorticity of the form q2iTS^/r = ±(2/a^) exp[-(r/s )^/a^)] at y/ s = 
±1 , z/s = 0, respectively. The coordinate r is measured radially 
from the vortex center and o is the vortex spread ( 0 = 0.5 for all 
computations) . The numerical mesh is specified for only y > 0 
since symmetries about y = 0 can be utilized. In addition to the 
uniform level of background turbulence, an isotropic Gaussian spot 
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of turbulence is positioned in the center of each vortex of strength 
(gp2irs/r)2 = 0.01 , such that <uu> = <w> = <ww> = (qg/3) exp[-(r/s) 2 / 
a 2 J] . A passive tracer of the same Gaussian character with strength 
C Q = 1 is also positioned over the vortex centers. The background 
constant density gradient [Fr 2rs 2 /r] 2 (g/po) 3p/3z = - 1 and the 
background axial velocity field is chosen such that the shear 
(2irs^/r) 3U/3z is that required to give the proper level of back- 
ground turbulence in accordance with super equilibrium theory. An 
upwash velocity is added to the solution in order to keep the vor- 
ticity in the center of the computational domain. This is equivalent 
to viewing the solution in a coordinate system that moves with the 
vortices. This upwash velocity is computed throughout the calcula- 
tion such that the centroid of the passive tracer is kept in the 
center of the domain. The initial conditions are shown schemati- 
cally in Fig. 4.2. 

The boundary conditions are given in terms of the moving co- 
ordinate system. The stream function is computed, from a Biot-Savart 
law as described in I with the addition of a contribution due to the 
moving coordinate system. The vorticity, density, and axial velocity 
are specified at the bottom of the domain by their values in the un- 
disturbed fluid; namely , vorticity is zero and the density and axial 
velocity are determined by the background gradients that were speci- 
fied initially. These three quantities are advected out through the 
top of the domain with the local vertical velocity as the domain 
descends into heavier fluid. At the outer edge of the domain, the 
condition on these quantities is that 3/3y = 0 . At the edges of 
the computational domain, except at the y = 0 boundary, turbulent 
quantities are assumed to have zero slope and the passive tracer 
is assumed to have zero value. The zero slope condition at the 
outer edge in the y direction is guite adequate at this boundary 
despite the presence of Brunt -Vaisaia waves since the computations 
were not carried so far in time that a reflecting wave became a 
problem. 

Nine cases were run with three different Froude numbers and 
three different values of initial background turbulence levels-. 

These runs are summarized in Table 1 below. 

All cases were run with a Reynolds number of 10^ and a fixed 
value of the integral scale A/s =0.1 . Results of cases 4, 5, and 
6 (no initial background turbulence) are shown in Figs. 4.3 and 4.4 
Figure 4.3 is a plot of the descent velocity of the vortex measured 
by tracking the centroid of the smoke, while Fig. 4.4 is a plot of 
the penetration depth of the vortices as a function of time. Simu- 
lations were carried out until the descent velocity dropped to less 
than 0.1 from an initial value of 0.5 . Note the much larger stop- 
ping times, as expected, for the higher Froude number (lower strati- 
fication) cases. It should be noted that the descent velocity is 
a monotonically decreasing function throughout the computation. 
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Figure 4.2. Schematic of initial conditions for the 
computation of a vortex pair descending 
in a stably stratified fluid. 




Figure 4.3. Descent velocity of the vortex pair versus time for cases 4, 5, and 6 
with Froude numbers 1, 5. and 10. 



TABLE 1. 



DESCENT OF 

A VORTEX PAIR INTO 

STABLE STRATIFICATION 

Run 

Fr 

(2tts/T) 2 q 2 

(2¥s/T)/(3U/3z) 

1 

1.0 

0.0090 

1.43 

2 

5.0 

0.0096 

0.50 

3 

10.0 

0.0091 

0.39 

4 

1.0 

0.0 

0.0 

5 

5.0 

o 

o 

0.0 

6 

10.0 

0.0 

0.0 

7 

1.0 

0.00098 

1.355 

8 

5.0 

0.00094 

0.310 

9 

1.0 

0.00091 

0.195 


For case 6, NY = 32 , NZ = 45 , = - 4.4 , z y = 4.4 , 

and yu = 7.0 . The grid in the z direction is uniform 
and the grid in the y direction is stretched to give 
maximum resolution near y = 0 . For all other cases; 

NY = 33 , NZ = 40 , z^ = - 3.0 , z u = 5.0 , and y u = 3.2 
with uniform grids in both directions. (Refer to Fig. 4.2 
for grid and domain size definitions.) 
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No acceleration of the pair is predicted. 

Results of cases 1, 2, 3, 7, 8, and 9 where background tur- 
bulence levels are specified are not shown since they are nearly 
identical to the respective Froude number cases. This observa- 
tion at first appears contradictory, but is simply explained by 
examining the maximum computed level of turbulence present in 
the vortex pair as a function of time. As can be seen in Fig. 
4.5, the internally produced turbulence in the descending oval 
far exceeds the ambient value and, therefore, it is this inter- 
nal process which intimately controls the descent of the pair 
in a stably stratified fluid. 

In order to measure the decay rate of the vortices, rolling 
moments induced on a follower aircraft of wing semispan sf have 
been computed and displayed in Fig. 4.6. It can be seen that the 
vortices have stopped their descent long before there has been 
any significant decay of rolling moment on the follower aircraft. 
This is a major finding in that vortices may, as a result of 
stable stratification, be observed to stop descending but may 
still be capable of inducing large rolling moments on an encount- 
ering aircraft. 


4.3 A Simple Estimate of the Penetration of a Vortex Pair into a 

Stably Stratif i e d Fluid, 

Analysis of the streamline pattern generated while under- 
taking the above computations has shown that the volume of fluid 
descending is "approximately" constant during much of the des- 
cent. This can be seen in Figs. 4.7 and 4.8 where streamline 
patterns are shown for the Fr = 1 and 10 case, respectively. 
During the final stages of descent, when the descent velocity 
is approximately 257. of the initial value, recirculating oval 
sizes do change but the additional penetration distance covered 
after this time is small when compared to the total penetration. 
This observation suggests a simple model which can be used to 
estimate maximum penetration. This model is based on an energy 
balance . 

The initial kinetic energy of the pair is related to the 
induced drag of the wing. Expressing the induced drag as a 
multiple of the drag of an elliptically loaded wing we write 

D ind. - (4 - 6 > 

where 8 is a function of the spread of vorticity a . For the 
sample computations above, 8 was numerically determined to be 
equal to 0.547 . The final potential energy of an oval of area 
A, carrying fluid with density equal to that at its initial 
height, is given by 
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Figure 4.5. Maximum internally produced turbulence levels 
as a function of tine for cases 4, 5, and 6 
(background turbulence level = 0) . 
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Figure 4.6. Rolling moment drop off as a function of time for follower aircraft 
of semispan s f . Cases 4, 5 and 6. 




Figure 4. 7 (b) . 


Streamline pattern at tr/2irs = 1 
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= 0.35 
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Figure 4.7 (c) . Streamline nattem at tr/2-rrs = 1.68 , 

contour interval = 0.06 , descent velocity 

dz 2-ts 
dt r 


0.08 . 
28 





Figure 4.8 (a). 


2 

Streamline pattern at tr/2irs = 0 , contour 
interval = 0.2 , descent velocity - ^ =0.5 . 
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Figure 4.8 (b) . Streamline pattern at 
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Figure 4.8 (c) . Streamline pattern at tr/2trs“ » 28 , 
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(4.7) 


P.E. 


max 



Ap in a linear stratification at the point of maximum descent 
is 


Ap = + 



max 


(4.8) 


With the assumption that the final oval area A is equal to the 
initial area, and that the oval area is given approximately by 
that of an ellipse having semi-axes of a translating point 
vortex pair 


A = tt ( 1 . 73) (2 . 09) s 2 (4.9) 

we find that the upper limit on the maximum descent of a vortex 
pair in a linear stratification, assuming constant oval area, is 


z 


max 

s 


1.65B 1/2 Fr 


(4.10) 


N 


is the Brunt-Vaisala 

N 


circular frequency in rad/ sec 




(4.11) 


Figure 4.9 shows a comparison between the maximum descent 
predicted by the theoretical upper limit and the results of a 
number of WAKE calculations. Both the approximate nature of 
this analysis and the fact that the numerical calculations con- 
tain self-induced turbulent dissipation, result in the discre- 
pancy between numerical computation and the energy penetration 
argument. This discrepancy increases as Froude number increases 
since the self-induced turbulence becomes a progressively more 
important factor in decelerating the pair the longer the des- 
cent period. The presence of ambient turbulence in the atmos- 
phere also places a limit on the maximum descent, independent 
of the stratification, see Donaldson and Bilanin Ref. 17. 

There, an approximate expression for the maximum descent of a 
vortex pair through a turbulent atmosphere is found to be 


z max _ 2 . 38 

s _ (q/V Q ) 

Thus, for an ambient turbulence level of 


2 ,..2 
9 /V, 


0.01 


(4.12) 

the maximum 
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Figure 4.9. Maximum penetration of a vortex pair in a linearly 
stratified incompressible atmosphere. 
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descent in an unstratified environment is approximately 
|S max /s l = 24 at Fr = » . 

The penetration model posed here appears to yield a simple 
conservative means of estimating the maximum descent of a vortex 
pair in a linearly stratified atmosphere. 
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5. VORTEX INTERACTIONS IN A WIND TUNNEL AND WAKE ROLL UP 


In this section two computations are described which demon- 
strate WAKE's ability to simulate flowfields in wind tunnels and to 
compute the roll up of a sheet of vorticity. 

5.1 Simulation of the Cross Flow in a Rectangular Wind Tunnel . 

It is anticipated that vortex wake flowfield data will be taken 
in wind tunnels with rectangular cross section v In addition, experi- 
ments have been carried out in NASA Langley's -Long Track to investi- 
gate vortex ground plane interaction. To demonstrate that the WAKE 
code can be adapted to compute these flowfield configurations, a 
sample calculation was undertaken of the ascent of two equal ' 
strength vortices of opposite sign in a square enclosure. 

The calculation simulates a square tunnel with walls located 
at y/s = ±1.0 , z/s = ±1.0 , and the flow is initialized with two 
Gaussian vortices of the form 

2-ns 2 n = ±4 0e -(r/s) 2 /(0.16) 2 (5.1) 

located at z/s =0 , y/s = ±0.31 , respectively, in a uniform 
background turbulence of magnitude (2irs/r)^q z = 0.01 . The turbu-, 
lent scale is initialized as A/s =0.2 , the Reynolds number is 10^ 
and the numerical grid is uniform in both directions with a spacing 
of Ay/s = Az/s = 1/16 . The law of the wall boundary condition, 
described in section 2, is applied at all four walls. Isopleths of 
the vorticity field at three times tr/2irs2 =0 , 5.0 , 9.0 are 
shown in Figs. 5.1 - 5.3. Note that the vortex pair is self-induced 
to move upwards and produces vorticity of opposite sign on the wind 
tunnel walls. Note also that this simulation could have been run 
using half the mesh points if the symmetry of the solution had been 
utilized. 

5.2 Roll Up of a Vortex Sheet 

The vorticity that is shed from the trailing edge of a wing is 
initially a vortex sheet which then rolls up into discrete vortices. 
The difficulties of computing sheet roll up by finite difference 
are tied to the problem of the high numerical resolution required 
to resolve the sheet of vorticity. To check to see if finite dif- 
ference is a feasible technique, a sample calculation has been per- 
formed. Assume there exists a parabolical load distribution on a 
wing such that 


r / r 0 = 1 - (y/s) 2 , | y | < s 


(5.2) 


35 




36 




.37 


Figure 5.2. Isopleths of vorticity, labels on the curves are ten times 
the fraction of the maximum value. 







Then the sheet strength is 


y = - dr/dy = t 


- 2yr o /s' 


0 


I y I < s 

otherwise 


(5.3) 


If it is assumed that the vorticity distribution is a Gaussian of 
width O) in the z direction, then 

/ h * o / / .. v 2 

■A 1 

2 1 
ns /r o = 


This vorticity distribution with (a ,/ s = 0.05 was used to initialize 
a laminar WAKE simulation with a Reynolds number of Re = 200 . The 
numerical domain was specified to be 0 < y/s < 2.0 , -2.0 < z/s < 
2.0 .with a resolution of Ay/s = 0.0625 and Az/s = 0.105 . The 
time : jsvolution of the vortex roll up is shown in Figs. 5.4 - 5.9 ' 

as isopleths of the vorticity in the cross flow plane. At a non- 
dimensional time of approximately 5, roll up is complete and an 
"axisymmetric" structure has evolved. Note the difference between 
the y/s and z/s scales. 


_£Z_ e~ ' z/( t' 


0 


I y I < s 

otherwise 


(5.4) 
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Figure 5.6. Vortex Sheet Roll Up - Isopleths of Vorticity 
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6. CONCLUSIONS AND RECOMMENDATIONS 


Several conclusions from this study can be made with regard 
to the ability of the atmosphere to dissipate aircraft vortex 
wakes . 

1. After vortex interactions in multiple vortex pair wakes 
have gone to completion, the resulting vortex pair in 
the absence of atmospheric effects is predicted to decay 
very slowly. Direct computation has shown that the cir- 
culation in a fixed contour and the rolling moment in- 
duced on an encountering aircraft decay asymptotically 
for large times as t “1/2 . 

2. The presence of ambient background turbulence in a 
neutral atmosphere has a strong influence on the decay 
rate of a vortex pair. Direct '.computation has shown 
that asymptotically for large time, the circulation in 
a fixed contour and the rolling moment induced on an 
encountering aircraft decay as t“2 . More detailed 
predictions of the decay rate must include specifica- 
tion of the integral scale of turbulence if comparisons 
are to be made with measured data. 

3. Crosswind shear results in the observed solitary vortex 
phenomenon. The vortex whose vorticity is of the sign 
opposite to that vorticity in the shear is dispersed. 

The decay of the solitary vortex is dependent on the 
amount of ambient turbulence associated with the shear. 
It is postulated that under stable conditions, where 
shear can be present without turbulence, the solitary 
vortex may persist for a long time. 

4. The classic picture of a vortex pair descending towards 
a ground plane following a hyperbolic trajectory is 
shown to be incorrect. Viscous interaction with the 
ground plane results in the generation and separation 
of secondary vorticity from the ground. This secondary 
vorticity reduces the strength of the vortex wake as 
well as results in a rebounding or bouncing of the pair 
away from the ground. 

5. A vortex pair descending into a stably stratified incom- 
pressible atmosphere is shown to have a finite penetra- 
tion dependent upon the severity of the stratification. 
It is suggested that proper accounting of turbulent 
transport is as important as proper accounting for the 
buoyancy produced vorticity in predicting the dissipa- 
tion of a vortex pair in a stably stratified fluid. A 
major finding is that once a vortex pair's descent is 
halted by stratification, the vortex is still capable 

of inducing a significant rolling moment on an encount- 
ering aircraft of a span less than or equal to the 
generator's span. 
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6. Lastly, sample calculations have shown that the WAKE code 
can be used to compute the initial sheet roll up as well 
as to simulate vortex flows in rectangular enclosures. 

This capability should prove useful in future wake 
related work. 

The results obtained from this study warrant the following 
recommendations . 

It is becoming increasingly apparent that the WAKE code has the 
potential of being an invaluable predictive tool in flow situations 
where it is impractical or even impossible to make measurements. 

The code's usefulness is, however, limited until proper validation 
against vortex related data can be completed. It is strongly 
recommended that a set of experiments be devised which will produce 
sufficient data appropriate to validate the WAKE code. 

Results of analytic studies reported here suggest that the 
atmosphere plays a powerful role in affecting the dissipation of 
aircraft vortex wakes. No comparison is directly made with data. 

The Department of Transportation, as part of their effort in the 
vortex hazard program, has monitored the descent of thousands of 
vortex wakes under various atmospheric conditions . It is suggested 
that a review of these data and correlation against predictions 
made herein will increase our understanding of the ability of the 
atmosphere to dissipate vortex wakes. 
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Abstract 

Enhanced dispersion of two-dimensional 
trailed vortex pairs within simplified 
neutral atmospheric backgrounds is studied 
numerically for three conditions : when the 
pair is imbedded in a constant turbulent 
bath (constant dissipation) ; when the pair 
is subjected to a mean cross-wind shear; 
and when the pair is near the ground. Tur- 
bulent transport is modeled using second- 
order closure turbulent transport theory. 
The turbulent background fields are con- 
structed using a superequilibrium approxi- 
mation. The computed results allow several 
general conclusions to be drawn with regard 
to the reduction in circulation of the vor- 
tex pair and the rolling moment induced on 
a following aircraft: (1) the rate of 

decay of a vortex pair increases with in- 
creasing background dissipation rate; (2) 
cross-wind shear disperses the vortex whose 
vorticity is opposite to the background; 
and (3) the proximity of a ground plane re- 
duces the hazard of the pair by scrubbing. 
The phenomenon of vortex bounce is explain- 
ed in terms of secondary vorticity produced 
at the ground plane. Qualitative compari- 
sons are made with available experimental 
data, and inferences of these results upon 
the persistence of aircraft trailing vor- 
tices are discussed. 

Nomenclature 

C passive tracer (C* = C/C o ) 

proportional to the induced rolling 
moment, Eq. (3) 

q root mean square of the turbulent 

kinetic energy (q* = q2rs/T o ) 

r radial coordinate measured from 

vortex center 

s initial semiseparation between 

vortices 

t time (t* = tr Q /2irs 2 ) 

U^ mean Cartesian velocity components 

<u i u-> ensemble averaged Reynolds stress 
J correlation 

Xj Cartesian coordinates 

z Q hydrodynamic roughness length 
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0 


6 normalized constant cross shear 

(2its 7 /r o )aV/3z 

r_ initial circulation 

o 

e turbulent dissipation rate 

(e* = (27r)3 s 4 £ / r 3) 

C streamwise component of vorticity 

(r = c27T S 2/r o ) 

A turbulent macroscale parameter 

v kinematic viscosity 

a vortex spread parameter 

ip streamfunction 

Introduction 

The decay of aircraft vortex wakes has 
been the focus of extensive investigation 
by both NASA and DOT, in an effort to de- 
termine hazard potential and wake lifetime. 
Because of the inherent difficulty in mak- 
ing measurements in the wakes of aircraft, 
investigations have included the develop- 
ment of computer codes to simulate vortex 
behavior and predict vortex interaction. 

The work described herein are results from 
one such code "VORTEX WAKE." 

Donaldson and Bilanin* present a re- 
view of the physics underlying vortex gen- 
eration and decay, and the approaches used 
to analyze these effects. A comprehensive 
review of the entire subject has recently 
been given by Hallock and Eberle. In 
Bilanin, Teske and Williamson* (hereafter 
referred to as "I"), we used the point- 
vortex approach of Rossow* to initiate a 
discussion of vortex merging, wherein vor- 
tices of like sign interact. The thrust in 
I was the application of turbulent trans- 
port theory to the merging process, to pro- 
vide for a complete description of the 
convective and dissipative processes within 
the vortices themselves. The two-dimen- 
sional, unsteady numerical code, "VORTEX 
WAKE," programs the Reynolds stress equa- 
tions modeled by Donaldson. 5 The formula- 
tion of the modeled equation and the 
numerical procedures involved in solving 
them are included in I, as is an extensive 
discussion of the merging of multiple 
trailing vortices and the minimum hazard 
potential for a wing. 

One important result from I is that a 
pair of oppositely-signed vortices exhibit 
a surprisingly slow rate of decay when 
descending in a calm, neutral atmosphere 
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away from the ground. These ideal condi- 
tions, fortunately, are not met often in 
practice. More commonly, vortex pairs en- 
counter background turbulence, mean wind 
gradients, or the ground while descending 
in an atmosphere which is not in general 
neutrally stable. The effects of back- 
ground turbulence, wind shear and ground 
interaction form the subject of this paper. 
Here, we will demonstrate through the use 
of our turbulent transport code that any of 
these quite naturally occurring effects may 
substantially reduce the lifetime of a vor- 
tex pair. Descent of a vortex pair in a 
stably stratified atmosphere will be the 
subject of a subsequent paper. 

In Section II, we briefly review the 
turbulent transport model and examine the 
superequilibrium limit of this model which 
is used to provide idealized atmospheric 
environments. Also, improvements in the 
numerical code are discussed here. 

In Section III, we examine the effects 
of background turbulence on a vortex pair, 
showing that decay increases with increas- 
ing dissipation rate. In Section IV, we 
investigate the effects of mean shear and 
offer an explanation of the ''solitary'' vor- 
tex phenomenon. Finally, in Section V, we 
investigate vortex interaction with a 
ground plane and explain the phenomenon of 
vortex bounce which has been observed on 
runways and in subscale tests. 


stress equations for 
written as : 


< U i Uj > 


may be 


production 

r iiK aiT' 

0 - - <u iV 5^ ' <U kV 3^ 

tendency towards isotropic 

isotropy dissipation 




( 1 ) 


where q^ = <u.u.> and A is the turbulent 
macroscale. Tfteproduction terms in (1) 
are balanced by a Rotta tendency-towards- 
isotropy term and an isotropic dissipation 
term. If the principle gradient is in the 
z-direction for the U velocity, for in- 
stance, then it may be shown that the solu- 
tion to (1) is 


q 2 

- 

2 (*m: 

i 2 


<uu> 

= 

(*S) : 

> 


<w> 

= 

<ww> = 

0.5 4 


<uw> 

= 

“0.354| 

(*!§: 

> 2 

<uv> 

= 

<vw> = 

0 



II . The Turbulence Model and "Vortex Wake" 

^ ot j e 

A need to predict turbulent transport 
in rapidly rotating, as well as strati- 
fied flows, led Donaldson to develop a 
second-order closure model of turbulent 
transport. A second-order model implies 
that partial differential equations are 
solved for the Reynolds stress tensor 
<u^uj> as well as for the mean ensemble 
averaged velocity components. Higher- 

order unknowns and certain second-order 
correlations in the Reynolds stress equa- 
tions are replaced by modeled terms con- 
taining constants and second-order correla- 
tions. These constants are evaluated by 
comparison with measurements made of funda- 
mental fluid flows. A complete summariza- 
tion of the current state of incompressible 
modeling as used here may be found in 
Lewellen and Teske.® The computations 
which follow idealize the neutral atmos- 
phere by assuming constant shear and homo- 
geneous turbulence in the absence of a vor- 
tex wake. The background turbulence levels 
are determined by the "superequilibrium" 
limit of the turbulent transport model. ^ 
Superequilibrium limit is that turbulent 
fluid state which is approached at high 
Reynolds number when time rates of change 
and diffusion of second-order correlations 
are negligible. When the diffusion and 
rate terms are set equal to zero in the 
rate equations for the second-order velo- 
city correlations and a high Reynolds num- 
ber is assumed, the resulting Reynolds 


In the vicinity of a ground plane , the 
velocity gradient is inversely proportional 
to z , (with z measured upwards from the' 
ground) while the macroscale is proportional 
to z , maintaining a constant product of 
A3U/8Z . 

The mean and rate equations for the 
second-order equations correlations have 
been programmed into a two-dimensional un- 
steady incompressible computer code. The 
assumption of two-dimensionality implies 
that the phenomena to be investigated are 
such that axial gradients are insignifi- 
cant and can be neglected. Therefore, 
sinusoidal instability and vortex breakdown 
are phenomena which cannot be addressed 
with this code. 

The numerical solution technique (us- 
ing an ADI scheme) is given in I. Several 
extensions to the code since the publica- 
tion of I are the following. The mean 
equations have been reworked into stream- 
function-vorticity variables, eliminating 
the need to compute the pressure. A dir- 
ect solver developed by Swartztrauber and 
Sweet 7 determines the streamfunction, from 
which the velocity fields are obtained. 

In this way continuity is assured, even 
though a more difficult set of boundary 
conditions may sometimes result. To im- 
prove our simulation of the convective 
nature of the vortex flow problem, we now 
use a J6 representation of the convec- 
tive operator (Arakawa®) , wherein the con- 
servative and nonconservative operators are 
averaged. 
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For the specific computation of vortex 
pairs, an upwash velocity is added to keep 
the vortex pair in the computational domain. 
This is equivalent to sitting on the vortex 
pair as it descends. When the results are 
presented, however, the vertical coordin- 
ate is shifted to show the actual amount of 
descent which has resulted. The choice as 
to what defines the position of a vortex is 
arbitrary. In I the centroid of the vor- 
ticity distribution was used to define vor- 
tex location. This definition was adequate 
for vortex pair computations in the absence 
of a ground plane. As will be shown later, 
however, when vortices interact with a 
ground plane, vorticity of the opposite 
sign is generated at the ground. This ad- 
ditional vorticity limits the 'usefulness of 
the centroid concept to fix position. To 
overcome this difficulty a passive tracer 
(neutral smoke) is initially prescribed in 
the center of the vortex and the maximum 
of this passive tracer is taken as the vor- 
tex center. As the pairs spread by turbu- 
lent diffusion, readjustment of the numeri- 
cal grid is undertaken to keep the 
boundaries away from the centers of vorti- 
city. This readjustment does not cause any 
numerical problems. 

As a convenience, initial distribu- 
tions of vorticity, second-order turbulence 
quantities and passive tracer are taken to 
be Gaussian. To reduce grid point require- 
ments, whenever possible, half-plane compu- 
tations are made. A uniform mesh is used 
except near the ground where additional ' 
resolution is, in general, needed to define 
the large shears that exist there. 

Ill . Vortex Pairs in Constant Turbulent 
Backgrounds 

In I the decay of a vortex pair in a 
constant dissipation rate turbulent bath 
was investigated to show that the dissipa- 
tion rate itself is not sufficient to de- 
termine the rate of decay of vortices. It 
was shown that in a neutrally stable atmos- 
phere, in the absence of mean cross wind 
shear, both the integral scale of the tur- 
bulent eddies, A , and the turbulent dis- 
sipation rate, e , must be specified to 
fix vortex wake decay rate. In this sec- 
tion, additional computations are described 
and simple decay laws are deduced from the 
computations . 

The computations are initialized by 
positioning a Gaussian of vorticity of the 
form t2irs^/r 0 = ± (2/o2) exp (-(r/sj’/a^) 
at y/s = ± 1 , z = 0 , respectively. 

The coordinate, r , is measured radially 
outward from the vortex center and o is 
the vortex spread, (c = 0.5 for all compu- 
tations) . The numerical mesh is specified 
only for y > 0 since symmetries about 
y = 0 can be utilized. In addition to a 
uniform background level of turbulence an 
isotropic Gaussian spot of turbulence is 
positioned in the center of each vortex of 
strength (q 0 2rs/r 0 )2 = 0.01 , such that 


2 2 

<uu> = <w> = <ww> = (q 0 /3) exp [-(r/s) 

/o^) . When a passive tracer is used, it 
begins- with the same Gaussian character and 
has strength C D = 1 . At the edges of the 
computation domain, except at the y = 0 
boundary, turbulent quantities are assumed 
to reach the superequilibrium state and 
tracer concentration is set equal to zero. 
The streamfunction is obtained from moment 
expansion of the Biot Savart law as des- 
cribed in I. The initial conditions are 
schematically shown in Fig. 1. 


z/s, W 



Figure 1. The initial conditions for the 
computation of the decay of a vortex pair 
in a neutrally stable turbulent atmosphere. 

As a baseline computation the ambient 
turbulent field was set equal to zero and 
the decay of a vortex pair in . a calm atmos- 
phere was simulated. In this computation, 
as in all computations presented here un- 
less otherwise noted, the integral scale 
parameter was initialized to be A/s =0.2 
and the Reynolds number was fixed at r D /v 
= 10^ . The evolution of the scale is 
controlled by a dynamic scale equation de- 
tailed in I. As can be seen by examining 
the contour plots of vorticity for several 
times as shown in Fig. 2, a vortex pair in 
a calm atmosphere is predicted to dissipate 
very slowly. The vorticity distribution 
retains its shape initially (except for 
some slight vortex-on-vortex straining) , 
and shows signs of trailing a wake after 20 
vortex flow times. The vorticity level has 
dropped a factor of 10 below its initial 
value after 100 vortex flow times. For a 
typical landing B-747 aircraft, a flow time 
roughly corresponds to about 10 seconds. 

The apparent persistence of the vortex pair 
is not in agreement with observation. 

A substantial increase in the decay 
rate occurs when the pair is placed in a 
turbulent flowfield. We assume that the 
turbulence is the result of a headwind 




Figure 3. Isopleths of vorticity for a vortex 
descending in a constant turbulent 1 
Notation explained in Fig. 2. 


gradient, 3U/3z , such that the back- 
ground turbulent level through (2) is non- 
zero and homogeneous . 

If the scale length is taken to be 
A/s = 0.2 and (3U/3z) / OV^s 2 ! = 0.56 
the dissipation rate, (2ir) 3 s**e/r3 = 0.0025 
is obtained. Isopleths of vorticity are 
shown in Fig. 3. In comparison with Fig. 2, 
it can be seen that substantial diffusion 
of the vorticity has occurred, with a drop 
in 'maximum value of an order-of-magnitude 
below the zero background case. In addi- 
tion, the structure of the vortex pair has 
been noticeably altered. It is clear from 
the shape of the isopleths at times greater 
than tT/2irs 2 = 20.0 that the vortex pair 
leaves a wake of low intensity vorticity. 

An extrapolation of the decay of the pair 
without background turbulence shows that 
over 10^ flow times must pass before the 
zero background pair would decay to the 
same level as this decayed pair in 100 flow 
times . 

Two quantities have been selected to 
quantify the amount of decay which occurs 
in a vortex pair; the circulation in a 
given area, and an integral which is 
approximately proportional to the rolling 
moment induced on a constant chord airfoil 
of semispan Sf . In Fig. 4 is shown the 



Figure 4. Circulation dropoff as a func- 
tion of time for square contours of size 
2e x 2e centered on the local maximum of 
vorticity for several background dissipa- 
tion ratio t* = (2u) 3s^e/rg . 

circulation computed about a square contour 
of dimension 2e * 2e centered on the maxi- 
mum of vorticity. The rate of decay as 
measured by r is independent of dissipa- 
tion rate for sufficiently large times.- 
These computations suggest that r - t” z , 
which is expected from simple diffusion 
arguments that treat vorticity as passive 


tracer. Curiously the e = 0 computation 
has a long time decay of T ~ t"^' 2 . This 
time behavior may be shown to be the diffu- 
sion limit of a turbulent decaying pair 
once mean fluid convection of vorticity 
is neglected. As a general expected obser- 
vation, when the integral scale is initial- 
ized to be the same between turbulent atmos- 
pheres of differing dissipation rate at 
comparable times after wake initialization, 
increased e results in decreased circula- 
tion computed over a fixed contour. By 
far, however; it is our inability to 
quantify A in the atmosphere which pre- 
cludes deriving a general decay curve for 
neutral turbulent atmospheres. 

Corresponding rolling moment predic- 
tions for a follower aircraft, having a 
coaxial encounter with the vortex, are 
given in Fig. 5. Here the rolling moment 



Figure 5. Rolling moment dropoff as a 
function of time for a follower aircraft of 
semispan sf/s = 1 having a coaxial vortex 
encounter within several background dissi- 
pation rates e* = (2u) 3 s 2, e/r3 . 

is proportional to 


■:/ 


/c +s f 


w (z c .y) (y - 


y c )dy 


(3) 


y c - s f 

and y and z are chosen so as to maxi- 
mize the integral at a given time. Typi- 
cally y c and z c are near the maximum 
In vorticity. The rolling moment curves 
for smaller follower aircraft are within 10 
percent of the curves shown at these e 
levels, and are thus fairly represented by 
the sf/s = 1 curves. Agajy^ the rolling 


moment curve decays as 
calm atmosphere and as 
lent atmosphere. 


for the 
t" 2 for the turbu- 


IV. A Vortex Pair Subjected to a Cross 
Flow Shear 
9 

Tombach has observed during aircraft 


*#. 
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Figure 6. Isopleths of vorticity of a vortex pair subjected to a cross shear. Values 
are ±5, ±2, ±1, ±.5 . 

Figure 7. Isopleths of turbulence for the pair in cross shear. Values .4, .3, .2, .1 
Figure 8. Isopleths of tracer for the pair in cross shear. Values .5, .2, .1, .05, .0 




flight tests that under certain atmospheric 
conditions the smoke in one vortex would 
vanish while the second vortex remained in- 
tact. This phenomenon has come to be cal- 
led the solitary vortex phenomenon. At the 
time of observation, it was conjectured 
that atmospheric shear might be responsible 
in some way. More recently, Rossow 10 using 
distributions of point vortices has shown 
that a vortex immersed in a cross shear is 
rapidly distorted when the shear is of op- 
posite sign of the vorticity. Computations 
were undertaken to see if the large distor- 
tion in vortex structure could generate 
turbulence and thereby diffuse the smoke 
seeded in the vortex to make it visible. 

The computation was initialized, as in 
Fig. .1, but now a cross shear of magnitude 
6 = (3V/3z) / (2rs 2 /r o ) = 0.5 was imposed. 
Figures 6 and 7 show isopleths of vorticity 
and turbulence. The vortex of sign oppo- 
site to the cross shear is distorted which 
results in an increased turbulent level in- 
side the vortex. Figure 8 is the computed 
dispersion of passive tracer which results. 
The fluid mechanisms which result in the 
enhanced dispersion of smoke in one vortex 
are now clear. It is concluded that it is 
cross wind that is responsible for the 
phenomenon of the solitary vortex. 

To quantify the decay of vortices in 
cross shear, the rolling moment was com- 
puted and is shown in Fig. 9. Note that 



rolling moments induced by the same vortex 
pair in a calm atmosphere. This general 
enhanced decay in the left vortex, is attrt 
buted to the turbulent diffusion associated 
with the cross shear. It has been specu- 
lated that in the atmosphere over periods 
of tens of minutes it is possible to have 
shear with turbulence when the atmosphere 
is stably stratified. One might speculate 
that under these conditions the solitary 
vortex may indeed persist for long periods 
of time. 

V. Vortex Pair in Ground Effect 

One of the most perplexing observa- 
tions of vortices has been the rebounding 
of a vortex as it approaches the ground. 

The phenomenon, we believe, was first 
observed by Dee and Nicholas H in 1968 and 
has been observed numerous times by resear- 
chers measuring vortex trajectories at 
airports. I- 2 In 1971 Harvey and Perry, 13 
prompted by the above observation, conduct- 
ed a series of tests in a wind tunnel with 
a moving floor to simulate correctly the 
ground plane. Using total head probes they 
concluded that the proximity of the vortex 
to the ground plane could result in boun- 
dary layer separation and secondary vorti- 
city which could explain the upward 
migration of the vortex pair. 

The computer code was modified to in- 
clude a viscous ground plane in an effort 
to simulate the ground interaction. The 
surface is fixed at a level plane z ■> z D 
(where z Q is the roughness height). As 
far as the streamfunction boundary condi- 
tions are concerned, this surface is merely 
another reflecting plane across which 
oppositively-signed image vortices influ- 
ence the flow (much as the image vortex of 
strength -r did in Sections III and IV). 
The vortex pair will be assumed sufficient- 
ly close to the surface so that the velo- 
city gradient, 3U/3z , maintaining a small 
turbulence level, may be taken to be 
inversely proportional to z . Since the 
scale length is proportional to z in the 
Monin-Obukhov sublayer, the background 
turbulence level may still be taken con- 
stant. All necessary turbulence quantities, 
however, will be given zero vertical gradi- 
ents at the hydrodynamic roughness height. 
We relate the vorticity level at the sur- 
face to the streamfunction there by a fin- 
ite difference representation of the 
Poisson equation and the velocity law-of- 
the-wall to give: 


Figure 9. Rolling moment dropoff as a 
.function of time for a follower aircraft 
of semispan sj/s = 1 having a coaxial en- 
counter. Normalized constant cross shear 
is 6 = 0.5. 
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the rolling moment induced by the right vor- 
tex is about half that induced by the left 
vortex at six flow times. As can be seen 
at six flow times, both the right and left 
vortices induce rolling moments on a fol- 
lower aircraft which are well below the 


where 1(4. is the streamfunction value at 
z+ = z r + iz , and z r is the reference z 
height above z 0 (typically z r /s = 0.01). 
When the flow simulation begins, (A) gives 
an immediate surface vorticity layer of op- 
posite sign to the vortex convecting toward 
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the surface from above. A sequence of con- 
tour plots of vorticity for the case where 
z Q /s = 10"^ , with the vortices initially 
positioned at y/s = ±1 ,z/s = 2 , is shown 
in Fig. 10. 



Figure 10. Isopleths of vorticity for a 
vortex pair descending toward a ground 
plane. Positive contours around the vortex 
follow from the given maximum values as 
given in Fig. 2. Negative contours gener- 
ated near the wall follow -5, -2, -1, 

-.05, -.02, -.01, -.005. 


As the vortex approaches the ground 
(its pair is across the reflecting plane at 
y = 0 ) , the boundary layer beneath grows 
and separates, resulting in a secondary 
vortex. The secondary vortex is of the 
opposite sign to the descending vortex 
around which it is advected. The second- 
ary vortex induces an upward and outward 
motion on the descending vortex pair, which 
results in a rebound (bounce) of the pair 
away from the ground plane. The rebound 
phenomenon described here suggests that the 
classical "inviscid” description of a pair 
of vortices following a hyperbolic traject- 
ory may not be a particularly good model. 

In Fig. 11 is shown the computed traject- 
ories of the maximum value of passive 
tracer . 



Figure 11. Trajectories of the maximum 
value of passive tracer carried by a vortex 
pair descending onto a ground plane. 


Note that a laminar computation has been 
undertaken at r D /v = 100 and the rebound- 
ing phenomenon still occurs. 

14 

Curiously, Barker and Crow have 
recently Investigated vortex pairs in 
ground effect and have attributed the re- 
bound to effects of finite vortex core 
radius. To demonstrate that this explana- 
tion is suspect, we have undertaken a com- 
putation with an inviscid boundary condi- 
tion applied on the ground (vanishing of 
vertical velocity). The vortex spread, 
o = 0 . 5 , is used and the results of the 
computation are summarized also in Fig. 11. 
As can be seen the rebounding phenomenon 
does not occur unless the viscous boundary 
condition is applied, even though the vor- 
tex cores are of finite size. It is seen 
that the finite difference solution with 
inviscid ground plane is in good agreement 
with the hyperbolic point vortex traject- 
ories. The slight drift upward results 
from diffusion of the smoke . 

The decay which results from scrubbing 
along a: viscous ground is presented in 
Fig. 12. Included is the rolling moment 
decay curve with zero background turbulence 
taken from Fig. 5. 
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Figure 12. Rolling moment dropoff as a 
function of time for follower aircraft 
having coaxial encounters with a vortex 
interacting with the ground. 

As can be seen, interaction with the ground 
provides a mechanism which greatly reduces 
the rolling moment induced on a follower 
aircraft over that induced by the same vor- 
tex pair left to decay out of ground effect. 

VI • Conclusions 

Two-dimensional aircraft trailing vor- 
tex pairs descending in a neutral atmos- 
phere are shown to have enhanced dispersion 
when the atmosphere is turbulent , when 
cross wind shears are present and when the 
vortices interact with the ground. These 
results are predictions from the finite 
difference computer code, "VORTEX WAKE , " 
which uses second-order closure turbulent 
modeling to include turbulent transport. 

In particular, the following conclusions 
are drawn with regard to neutral atmos- 
pheres: . 

1) The decay of a vortex pair immersed in 
a turbulent bath is enhanced with in- 
creasing dissipation rate providing 
the integral scale of the turbulence 
remains constant. In general, however, 
the dissipation rate, even in neutral 
atmospheres, is insufficient to 
quantify dispersing effectiveness of 
the atmosphere. 

2) The phenomenon of the solitary vortex 
is shown to result as a consequence of 
cross wind shear. The vortex having 
the sign of vorticity opposite that of 
the vorticity in the shear is convec- 
tively distorted, producing turbulence 
in the vortex. This turbulence then 
disperses the smoke which was intro- 
duced in the vortex for visualization 
purposes . 

3) The classic inviscid notion that a 
vortex pair descending towards a 
ground plane follows a trajectory that 


is hyperbolic is generally not correct. 
The vortex pair rebounds or bounces 
away from the ground as a result of 
separation of the ground boundary 
layer below the vortex. It is shown 
that this separation results in sec- 
ondary vortices which induce a velo- 
city on the descending vortices out- 
ward and upward, thereby explaining 
the rebounding phenomenon. It is 
speculated that when cross winds are 
present symmetries of the secondary 
vortices are destroyed. This results 
in tipping of the descending pair 
which has been observed. 
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